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Background 

Obvious need for cloud detection with satellite imagery 

• Improved real-time situational awareness 

• Mask to support product generation 

• Climate studies 

Robust VIIRS Cloud Mask (VCM) produced by NOAA 
with confidence indicators, phase, presence of aerosols, 
smoke, or shadows, and more 

• IR threshold tests, brightness temperature difference tests, 
solar reflectance tests, NIR and thermal IR thin cirrus tests 

• Inputs - Tbbs for M12, M13, M14, M15, M16, 14, 15; reflectance 
for Ml, M4, M5, M7, M9, M10, Mil, II and 12 

• Variety of dynamic ancillary data -surface T, PW and wind 
speed, snow/ice cover - hundreds of adjustable parameters 

• Heritage - Saunders and Kriebel (1988), Stowe et al (1995) for 
AVHRR and Ackerman et al. (1997) for MODIS, and Gustafson et 
al. (1994), Hutchison et al. 2005 (for oceans) 
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Motivation 


VIIRS cloud mask already exists So why another? 

VCM design to meet NOAA needs, may not address 
broader community needs like: 

• Easy to implement approach for application real-time and 
case study 

• Tuning by user for better regional and local performance 

• Cross platform (instrument) consistency of approach) 

Two channel (11 pm and 3.7 or 3.9 pm) / three test 
with dynamic spatially and temporally varying 
thresholds may provide flexibility to the community 

Approach adapted from that of Haines et al. 2005; 
Jedlovec et al., 2008; Jedlovec 2009 
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Drivers to the Approach with VIIRS 

Easy to apply, accurate in variety of situations 

Day and night applications for diagnostic and climate 
studies, tunability for local applications 

Algorithm consistency between satellites 

* Rely on simple infrared threshold test in 11 pm window 
channel and emissivity differences between surface and 
various cloud types at 11 and 3 pm as represented in 
channel brightness temperatures 

Tbll - cold for high, thick clouds 

* Emissivity clouds is lower at 3 pm than 11 pm providing 
lower brightness temperatures at night but potentially 
larger (during day) when reflected solar radiation is present 

Tbll - Tb3 » 0 clouds at night 

Tbll - Tb3 « 0 clouds during day (solar illumination) 
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Cloud Tests 


Split window spectral test 

• If {Tbll(x,y,t)-Tb3(x,y,t) < 0 and < [(CND(x,y,t')-AND]} ^ CLOUD 

• If {Tbll(x,y,t)-Tb3(x,y,t) > 0 and < [(CPD(x,y,t')+APD]} ^ CLOUD 

AND - negative difference adjustment factor (4.0) 

APD - positive difference adjustment factor (2.5) 

Longwave infrared channel test 

• If {Tbll(x,y,t)-CWLW(x,y,t') < -ALW} ^ CLOUD 

ALW - longwave IR adjustment factor (18.5) 

Adjustment factors are subjectively determined to tune 
the algorithm performance for regional applications 
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Composite Thresholds 


Appropriate thresholds often change with season, 
time day, and geographic region 

Use imagery over previous 7-20 to determine 
appropriate threshold on a pixel by pixel basis 

• Largest negative and smallest positive Tbll - Tb3 pm 
difference image over the previous 20 days assumed to 
provide estimate of clear sky difference value - difference 
image values smaller than these thresholds probable 
correspond to cloud-free conditions 

• Warmest Tbll over previous 20 days at corresponding 
observation time represents cloud-free conditions 

Resulting composite images form spatially and 
temporally (recalculated each day) varying 
thresholds for cloud detection 
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Composites Use in Tests (Night) 
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Satellite Imagery / Cloud Mask (Day) 
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Satellite Imagery / Cloud Mask (Day) 
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Satellite Imagery / Cloud Mask (Night) 
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Satellite Imagery / Cloud Mask (Night) 



By 



1 * -f .' - 




- — 


, rvwii 

I W.’, 

i % S3 

#‘ 


■ 


A 

\ 
■ ~s 

. -£ f- 



12/15/2014 


IN12A-04 -- A simplified Approach to Cloud 
Masking with VIIRS in the S-NPP / JPSS Era 


12 




Verification / Validation 


How good / accurate is the cloud mask? 

Compare with other cloud masks, other satellite 
sensors (e.g., Calipso) and / or against ground truth 

Ground truth is difficult to obtain 


Develop ground truth data set for specific regions with 
trained remote sensing scientists / students 

• 5 weeks (October-Nov) , day / night passes over CONUS, total 
number of images and points 

• 3-5 satellite remote sensing experts manually look at imagery 
to determine % cloud cover, 30x30km (pixel) region - used 
visible and infrared imagery 
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Cloud Mask Validation Points 

CONUS satellite image with validation points 
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Cloud Mask Validation Points 


CONUS satellite image with validation points 
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Statistics (VI IRS) 


Performance against 42 validation points over from 
October 1 - November 5, 2014 



OBS 

C 

1 

F 

U 

Z 

H 

NUMBER 

3893 

2892 

1001 

771 

230 

993 

1899 

PERCENT 

100% 

74% 

26% 

20% 

6% 

26% 

48% 


C (correct) - total points and percent correct 
I (incorrect) - total points and percent incorrect 

F (false alarms) - the clear points identified as cloudy or the over-determination of clouds 
U (misses) - undetected cloud points or an under-determination of clouds 
Z (correct clear) - the points correctly determined to be cloud-free 
H (hits) - the points correctly determined to be cloudy 




VIIRS 

HR (hit rate) 

= H / (H + U) 

89.1% 

FAR (false alarm rate) 

= F / (F + Z) 

43.7% 

SS (skill score) 

= HR - FAR 

0.455 
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Performance 
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Spatially and temporally varying thresholds 
provide situational ... 

Tuning parameters - land water, etc. 
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Summary / Next Steps 


• Verification results for 5 weeks in Fall 2014 indicate an overall 
hit rate of 89.1%, however, the false alarm rate is quite high at 
43.1% suggesting a substantial over-determination of clouds 

• Preliminary results suggest that "adjustment factors" need 
tuning to meet past algorithm performance with MODIS 
imagery 

Next steps 

• Tune adjustment parameters 

• Validate in multiple seasons 

• Compare to NOAA VCM with external validation 

• Extend beyond CONUS 

• Global verification 
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Backup Charts 
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